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The competition between magnetism and Kondo effect is the main effect determining the phase 
diagram of heavy fermion systems, ft gives rise to a quantum critical point which governs the 
low temperature properties of these materials. However, experimental results made it clear that a 
fundamental ingredient is missing in this description, namely superconductivity. In this paper we 
make a step forward in the direction of incorporating superconductivity and study the mutual effects 
of this phase and antiferromagnetism in the phase diagram of heavy fermion metals. Our approach is 
based on a Ginzburg-Landau theory describing superconductivity and antiferromagnetism in a metal 
with quantum corrections taken into account through an effective potential. The proximity of an 
antiferromagnetic instability extends the region of superconductivity in the phase diagram and drives 
this transition into a first order one. On the other hand superconducting quantum fluctuations near 
a metallic antiferromagnetic quantum critical point gives rise to a first order transition from a low 
moment to a high moment state in the antiferromagnet. Antiferromagnetism and superconductivity 
may both collapse at a quantum bicritical point whose properties we calculate. 



I. INTRODUCTION 

Competition between superconductivity (SC) and 
magnetism plays an important role in determining the 
physical properties of strongly correlated superconduct- 
ing materials, such as high-T c superconductors^ and 
heavy fermions^^iSii. The study of the phase dia- 
grams of such materials suggests that these phases are 
intrinsically related. This relation and the interface be- 
tween these phases are among the fundamental issues 
that have not yet been clarified in condensed matter 
physics. In particular, the case of high-T c oxides rep- 
resents a formidable problem as the nature of the normal 
phase itself is not understood. 

The study of heavy fermion materials has focused until 
recently on the competition between long range magnetic 
ordering due to the RKKY interaction and the Kondo 
effect. However, low temperature experiments in these 
systems have shown that they can exhibit superconduc- 
tivity near or in coexistence with an antiferromagnetic 
(AF) phase close to a quantum critical point (QCP) 8 . 
This can be conveniently tuned by varying some element 
concentration or pressure^. In distinction with the high- 
Tc materials in the case of heavy fermions we have a quite 
reasonable understanding of their normal phase. Below 
the coherence line T co h this is essentially a strongly cor- 
related, nearly antiferromagnetic, Fermi liquid^. Above 
this line and in particular at the QCP, the non-Fermi 
liquid behavior is well understood in terms of quantum 
criticality and the associated critical exponents^. 

In this paper our aim is to investigate the mutual 
effects of superconductivity and antiferromagnetism at 
zero temperature and how this determines the phase di- 



agram of three dimensional heavy fermions (d=3). Our 
knowledge of the metallic phase in these materials pro- 
vides a solid starting point for the approach we use. We 
consider a Ginzburg-Landau functional which contains 
the magnetic and superconducting order parameters and 
their coupling. We apply the effective potential method 
used in quantum field theorjiiii^ to calculate the quan- 
tum corrections to this classical action. This method 
is implemented up to one loop order and requires the 
knowledge of the propagator solutions of the quadratic 
functionals associated with the decoupled fields. In such 
approach, the normal, metallic paramagnetic phase is de- 
scribed by the usual dissipative propagator that takes 
into account the dynamics of the spin fluctuations (para- 
magnons) in the nearly antiferromagnetic metaU2*i4. 

We investigate several phase diagrams in d = 3. A pos- 
sible scenario is that of a quantum bicritical point sepa- 
rating a metallic antiferromagnet from a superconducting 
phase whose universality class has been identified^. The 
model also describes the case where the AF and SC tran- 
sitions occur at distinct quantum critical points (QCPs), 
so that, there is a normal state between the AF and SC 
phases. The proximity to an AF instability drives the 
SC transition to the normal state into a first order one 
and enlarges the region of the phase diagram where su- 
perconductivity exists^. On the other hand supercon- 
ducting fluctuations drive an AF-QCP into a discontin- 
uous transition and give rise to a new AF phase with a 
reduced value of the order parameter. We present the 
Ginzburg-Landau model and the calculation of the ef- 
fective potential in detail since, as far as we know, this 
method has not been sufficiently explored in condensed 
matter physics. Although the model refers specifically to 
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the coupling between AF and SC order parameters, the 
results obtained should be considered more general, inde- 
pendent of the specific nature of the instabilities, as long 
as, the dynamics are similar. For example, a structural 
instability associated with a soft optical phonon mode 
has a dynamic description which is similar to that given 
to the superconductor here. 

It is interesting that recent experiments in heavy 
fermions suggest the existence of intrinsic inhomo- 
geneities in these materials near their magnetic quantum 
critical points. Our model provides a natural explana- 
tion for these inhomogeneities, attributing them to the 
weak first order transitions and the associated phenom- 
ena of phase coexistence and spinodals arising from the 
quantum corrections due to the coupling between differ- 
ent order parameters. 

We point out that similar models describing the com- 
petition between superconductivity and antiferromag- 
netism have been proposedi£*i2*i*i. The emphasis how- 
ever has been on the classical aspects of this phenomenon 
and their aim was to describe high-T c superconductors. 
The validity of this approach for this class of systems is 
now subject of intense debat o 20 ' 21 . 

II. PARAMAGNON AND 
SUPERCONDUCTING PROPAGATORS 

The model Ginzburg-Landau action contains three real 
fields. Two fields, fa and fa, correspond to the two 
components of the superconductor order parameter (the 
ground state wave function, for example). The other 
field fa, for simplicity represents a one component an- 
tiferromagnetic order parameter. The results are imme- 
diately generalized to a three component AF vector field, 
with the unique consequence of changing some numerical 
factors as discussed below. In order to include quan- 
tum fluctuations by the effective potential method it is 
convenient to describe the quadratic parts of the action 
by the associated propagators. The propagator associ- 
ated with the free, quadratic action of the superconduc- 
tor can be directly obtained from a quantum Ginzburg- 
Landau action 22 . However, finding the appropriate form 
of the superconductor propagator is a difficult problem 
since the nature of the fluctuations driving the supercon- 
ducting transition can change considerably the quadratic 
part of this generalized action. For metallic host where 



pair breaking interactions due to magnetic impurities de- 
stroy superconductivity the quadratic action has been 
obtained in Refi^ and is associated with a dynamic ex- 
ponent z = 2. The case of BCS superconductors in which 
any attractive interaction U makes the system supercon- 
ducting at T=0 has been studied in Ref. 24 . In this case 
the form of the action directly reflects the fact that the 
quantum critical point occurs for {7 = 0. Here we are in- 
terested in pure systems, so that criticality is achieved by 
pressure and impurities are not considered. Also, in our 
case, superconductivity is certainly non-BCS and finding 
the correct propagator is in itself a very difficult prob- 
lem. In order to make progress, we choose to consider the 
simplest generalization of the classical Ginzburg-Landau 
action, 

G (k) = G (cj,q) = 1 (1) 

where A; is a four-vector (u>, q) with k 2 = oj 2 — q 2 and 
we took the Fermi wave vector qp = 0. This propa- 
gator assumes the existence of a gap or pseudogap in 
the excitation spectrum of the phase precursive to the 
superconducting one. There is plenty of evidence for 
such pseudogaps in High-Tc superconductors 2 ^ and also 
now in heavy-fermions systems as shown in recent ex- 
periments^. Attractive Hubbard models also show nor- 
mal phases with charge ordering and a gap for excita- 
tions that vanishes in the superconductor quantum crit- 
ical points. 

We work in Minkowsky space since the present super- 
conductor propagator is Lorentz invariant and we want to 
obtain the effective potential as in quantum field theory, 
following closely the work of Coleman and Weinberg 12 . 
It's clear from Eq. that time acts as an extra direc- 
tion on the same footing as the spatial ones, so that, the 
associated dynamic exponent is z = 1. 

The quadratic functional associated with the magnetic 
part represented by the field fa, the sub-lattice magne- 
tization, takes into account the dissipative nature of the 
paramagnons near the magnetic phase transition. This 
approximation to the full paramagnon propagator can 
be found in detail in the work by Herta 1 ^ and here we 
present a brief derivation. The expansion of the mag- 
netic free-energy functional in a power series of the order 
parameter can be written as 



S mag (fa) = / <fkv 2 {k)\fa{k)\ 2 + 



d kid k 2 d k 3 d k i V4(ki,k 2 ,k 3 ,k i )fa(ki)fa(k 2 )fa(k 3 )fa(k4,)5(k 1 +k2 + k 3 + k i ) + 



(2) 



r 



and the coefficients v m are proportional to a loop of m electron propagators 2 ^. For a Hubbard interaction U the 
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susceptibility \ calculated in RPA is 



Imx(w, q) = Im 



1 - Uxo(u,q) 



(3) 



where xo is the susceptibility of the free-electron model 29 . 
This susceptibility gives a good representation of the 
paramagnon propagator associated with the quadratic 
coefficient i>2- If we are near a magnetic instability we 
can expand xo for long wavelengths and low frequencies. 
Substituting this expansion in the RPA form of the sus- 
ceptibility Eq. J3J and choosing the units appropriately 
we can write the AF paramagnon propagator as 



D ((J,q) 



i\u}\r - 



(4) 



where r is a characteristic relaxation time and m 2 is re- 
lated to local Coulomb repulsion U and the density of 
states at the Fermi level N(Ep) by 



m 2 = 1 - U N(Ep). 



(5) 



The quadratic parts of the effective functional are then 
completely characterized by the propagators of Eq. (|T|) 
and Eq. (J3J. Notice that the magnetic propagator, 
Eq. Q), is sufficiently general to be associated with a 
Hamiltonian which preserves rotational invariancei 3 -. In 
this case 03 is replaced by a vector and in the second 
term of Eq. the quartic interaction is written in the 
form (|</> 3 1 2 \<p 3 1 2 ) as appropriate for vector fields. The 
use of a vector field modifies the effective potential only 
by numerical factors as we show in the next sections. 



III. CLASSICAL POTENTIAL 



The classical part of the potential is given by, 

V cl 



, 3 ) = -m 2 (0? + ^) + -m^+ 



<h) + V p (<h) + Vi(4> u <h,<h) (6) 



where the self-interaction of the superconductor field is 

A 



V, 



and that of the antiferromagnet is given by, 



= 4^3- 



(7) 



(8) 



The coupling g above is related to the coefficient U4 in 
the expansion of Eq. with all dependence on the four- 
vector k ignored. Finally, the last term is the (minimum) 
interaction between the relevant fields, 

Vi(th,<h,<h)=u(<ft + <&)<&■ (9) 

This term is the first allowed by symmetry on a series ex- 
pansion of the interaction. It is possible that in regions 



of the phase diagram where phase coexistence between 
superconductivity and antifcrromagnetism occurs higher 
order terms should be included in the free energy. How- 
ever, we will not consider the case of coexistence in the 
present work. Notice that for u > 0, which is the case 
here, superconductivity and antiferromagnetism are in 
competition and the possibility of AF acting as a pair 
breaking perturbation is considered. 

The minima of the classical potential yield the classical 
ground states of the systemic. The different possibilities 
for the phase diagrams described by the classical action 
at mean field levelE'^ 1 ?.?? are shown schematically in 
Fig. ^ The T = transitions are all continuous, sec- 
ond order phase transitions. Notice that in Eq. © there 
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FIG. 1: Possible classical temperature phase diagrams for a 
heavy fermion system (schematic). AF and SC refer to the 
superconductor and antiferromagnetic phases respectively. 



are two mass terms, m 2 and m 2 , which represent respec- 
tively, the distance to the superconductor and the anti- 
ferromagnetic quantum critical points. These two terms 
appear in the quadratic parts of the action and are al- 
ready included in the propagators of Eq. (JIJ and Eq. Q 
to be used in the computation of the effective potential 
below. 

The classical potential can be generalized to a rota- 
tional invariant form replacing the scalar field 4>3 by a 
vector. As discussed before the self-interaction is then of 
the form 



V v 



|l0sl 2 



103 



i,3 



(10) 



where i and j label the components of the vector. For 
this generalized potential, however, the minima depend 
only on the modulus of the vector </> 3 and the use of 
this vector field brings only slight modifications to the 
effective potential. 

In the next sections we obtain quantum corrections 
to the classical potential. These corrections can modify 
significantly the physics of the problem and the phase 
diagrams changing, for example, the order of the transi- 
tions. 



IV. ONE LOOP EFFECTIVE POTENTIAL 

The first quantum correction to the potential can 
be obtained by the summation of all one loop dia- 
grams (Fig.0). 



4 




FIG. 2: One loop diagrams. The superconductor fields are 
represented by a or (3 = 1,2. The doted line represents the 
unusual propagator of Eq. @. 



We apply the general method proposed by ColemanSi 
with minimum modifications to account for the different 
nature of the propagators in our problem. The sum over 
the field indices can be easily done if we define a vertex 
matrix M , given by 



[M] 



lm 



°50 ; d0„ 



M={<M 



11) 



and then take the trace. In Eq. (|llfl the propagator 
(K l = G or D ) of Eq. {IJ or Eq. (gj is incorporated 
in the definition of the matrix. We draw the loops with 
arrows and choose the outgoing propagator of each ver- 



tex to be included in the associated element. The matrix 
M is then obtained deriving the classical potential with 
respect to the fields {0} and taking the values of these 
derivatives at the classical values of the fields, {0; c }. The 
sum of diagrams with the correct Wick factors is formally 
done in momentum space and using the property of the 
trace 



Tr[ln(l - M)] = lndet[(l - M% 



we get 



h / d 4 /clndet[l-M(fc)] 



(12) 



(13) 



The 3x3 matrix M can be simplified if we choose the 
classical minimum of the superconductor fields imposing 
4>2c = (this can be done because the minimum depends 



k 2 : 



only on the modulus 
clidean space, so that, K" = to" 
the first quantum correction can be written as 



J. Hence, rotating to Eu- 
+ q 2 and using h = 1 units 



^ (1) (01c03c) = ~ 



d 4 k 



In 

D 2 



k 2 

he) 



(k 2 + m 2 )(\uj\T + q 2 + in 2 ) 



In 



B( 



k 2 



1 C((j) lc ,(j) 3c ) \ 



NT 



(14) 



where 

A(& CJ &c) = (A/6)0? c + 2u0§ c (15) 

S(^ic,03c) = (A/2)0f c + 2nC (16) 

C(cl> lc ,cf> 3c )=2u4>l c + {g/2)<Pl c (17) 

£>(0ic,03c) = 4u0i c 3 c (18) 

The total effective potential with first order quantum cor- 
rections is then given by 

K/(01c, 03c) = Wlc, 03c) + V {1) (01c, 3c ) (19) 

where V c i is the classical potential of Eq. © and is 
the first quantum correction of order H of Eq. (|14l) . In the 
subsequent sections we consider the effects of these quan- 
tum correction for the classical results. We give special 
attention to the possibility of fluctuation induced symme- 
try breaking and first order phase transitions at T = 0. 
We point out again that the results for an isotropic spin 
model are similar since, when dealing with a vector mag- 
netic order parameter </> 3 , we can use the fact that the 
classical minima depend only on the modulus of this vec- 
tor. The modifications which arise in the final effective 
potential are simply different numerical factors coming 
from new but equivalent terms in Eq. I|14|) . 



V. POSSIBLE PHASE DIAGRAMS 

The ground state phase diagram of a heavy-fermion 
with AF and SC phases can be obtained from the classi- 
cal action by varying the masses m and m p . The normal 
paramagnetic state has m 2 > and m 2 > and for these 
values of the masses the configuration 0i c = 03 C = min- 
imizes the action. The AF and SC states have m 2 < 
(0 3c ^ 0) and to 2 < (0 lc ^ 0), respectively. The T = 
transitions between the different phases, which are tuned 
varying the masses, are up to the classical level all con- 
tinuous, second order transitions. The different phase 
diagrams, at mean field levelSiiLSI were shown already 
in Fig. n However, renormalization group studies of the 
classical action yield that a bicritical point, as that of 
Fig. ^3 but occurring at finite temperatures, is stable 
only if the number of components of the superconductor 
(ns) and antiferromagnet (uaf) order parameters^ are 
such that n = us+uaf < 4. For n > 4 the SO(n) bicriti- 
cal point at T ^ becomes unstable and a phase diagram 
like that of Fig.Qp may occur. The above constraint im- 
plies, for example, that in a heavy fermion antiferromag- 
net of the Ising type, long range magnetic order can not 
coexist with superconductivity (0i ^ and 03 ^ 0). No- 
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tice that quantum corrections may destroy the classical 
SO(n) symmetry associated with the quantum bicritical 
point in Fig. due to the different dynamics of the 
magnetic and superconductor fluctuations. 

In the next sections, we study the quantum effects in 
two of the possible phase diagrams of Fig. ^ namely, 
cases (A) and (B). We begin considering case (A) where 
there is a normal phase separating the SC and AF phases 
and the transitions occur at different QCPs. The quan- 
tum effects become important in the normal region near 
both SC and AF phases where symmetry breaking and 
fluctuation induced quantum first order transitions can 
occur. We distinguish here between the case the system 
becomes antiferromagnetic in the presence of supercon- 
ducting fluctuations from that where it becomes super- 
conductor in the presence of antiferromagnetic fluctua- 
tions. 

Case (C) of Fig. ^ with some additional assumptions, 
can also be described in our approach. In this case 03 
should be identified as the longitudinal component of the 
sub-lattice magnetization with a relaxational dynamics 
described by the propagator, Eq. Also we have to 
neglect spin wave excitations which, at T = 0, could be 
imagined as quenched if, for example, there is a suffi- 
ciently large anisotropy gap in their spectrum. In any 
case this situation has revealed up to now analytically 
intractable. Although it is possible to perform an ex- 
pansion of the last logarithmic term in Eq. I|14l) in pow- 
ers of u, this would be incompatible with the one- loop 
expansion which sums all powers of this quantity when 
performing the relevant integrals analytically, as we do 
below. 



VI. QUANTUM EFFECTS IN THE NORMAL 
PHASE BETWEEN THE AF AND SC-QCPS 



state, we look for a partially symmetry broken phase with 
0i c ^ but 03c = 0. In this case D((f>i c ,(j) 3c ) = and 
the quantum correction given by Eq. (|14|) can be written 
as 



+ 



d 4 k 

d A k 



In 1 



2 J (2tt) 4 
d 4 k 



In 1 



In 1 



k 2 + to 2 

k 2 + to 2 
2ucj)\ c 



(20) 



The first two integrations depend only on the modulus 
of the four-dimensional vector k and time enters as an 
extra dimension as we are dealing with a Lorentz in- 
variant case. This arises since the QCP of the super- 
conductor transition (SQCP) has an associated dynamic 
exponent z = 1. Therefore a cut-off regularization can 
be done as usuaU^. However, in the last integration we 
have anisotropy between time and space since the dy- 
namic exponent which characterizes the scaling of time 
takes the value z = 2. Hence, if we use a cut-off A to the 
momentum, the correspondent frequency cut-offl^ must 
be A z = A 2 . The integrations in Eq. (|20|l can be easily 
performed and the results expanded in powers of the pa- 
rameter m 2 supposing proximity to the superconductor 
transition (m 2 ~ 0). The renormalization is done begin- 
ning from the massless case m 2 = and generalized for 
small m 2 following closely the renormalization procedure 
for the charged superfluid 2 ^. The effective potential 



Fluctuation effects in the normal phase near 
superconductivity 



, 03c = 0) = V cl (cf> lc , 03c = 0) + F« (0 lc ) (21) 



We are interested here in the superconducting-normal 
quantum phase transition and therefore, near the SC 



where V c i is the classical potential given by Eq. © , is 



J 



V e 



-TO 0j 



4! 



(2tt)< 



15 



3 j. 2 
;TO p U0 lc 



r 



(22) 



Terms proportional to A 2 and to 2 A were neglected since potential can be written as a function of its extremum, 

they must be much smaller than the classical term pro- (0), which determines if the system is in the normal or 

portional to A in the small coupling limit. The effective in the broken symmetry superconducting phase. Higher 

potential, Eq. (|22|l . contains only the lowest order term in powers in u are not neglected since we have no classical 

to 2 since we are near both transitions and consequently term proportional to u. Consequently, even if m 2 > 0, 

to p is also small. The term between brackets is the first balancing the classical A term with the u terms of the 

quantum correction of order h. As usual, the effective quantum correction, it is possible to obtain asymmetric 



6 



minima for <fix c , i.e., a symmetry breaking in the normal 
state induced by the coupling u. For small masses m p , 
independently of the value of u, the asymmetric minima 
are close to the origin (\>\ c — 0, as we can show by nu- 
merical inspection of the full effective potential (this also 
justifies the expansion for small m 2 ). Furthermore, with 
all minima close to the origin, we can neglect </> 5 terms 
that cause instabilities to the potential (fortunately these 
terms become relevant far from the origin and hence away 
from the region of validity of the loop expansion) . 

For convenience, in Eq. Ill'2|l we have introduced the 
separation A 2 = m 2 + m 2 > (see Fig. between 
the magnetic and superconductor T — transitions and 
study the phase diagram as A 2 is reduced. Notice that 
A 2 measures the distance between the second order mean 
field QCPs in Fig. ^ and is also taken as a small quan- 
tity besides m 2 and m 2 . Now, studying the minima of 
the effective potential we see that the quantum correc- 
tions induce symmetry breaking and the region where 
superconductivity is found in the phase diagram is ex- 
tended. The shift of the superconductor QCP occurs to- 
wards the antiferromagnetic quantum critical point (AF- 
QCP), but in the paramagnetic phase. The analysis of 
the extrema of the potential can be easily carried out 
and has been presented before^. We also find that the 
quantum fluctuations changes the nature of the supercon- 
ducting transition to a first order one. The latent heat 
and spinodal points for this transition have also been cal- 
culated^. Between the two spinodals there is an interval 
of coexistence of the superconducting phase with regions 
of strong antiferromagnetic fluctuations due to the oc- 
currence of metastablc minima. 



B. Fluctuation effects in the normal phase near the 
AF-QCP 

We follow the same procedure here as in the last sec- 
tion. Now we are interested in the AF quantum phase 
transition and therefore we look for a partially symme- 
try broken phase with 4>^ c ^ but <\>\ c = 0. We have 
once again D(tf>i c , (f>^ c ) = in the quantum correction 
given by Eq. (|14|l and the three resulting integrals have 
the same functional form of Eq. i|20|l but with different 
4> c dependence. Integration and renormalization can be 
done as before beginning from the m 2 = case. The full 
effective potential 



0. 



) = V C 



l{<Plc 



0,<hc) + VW(<hc) (23) 



is similar to that 
cept for the renori 



be neglected since 
sical potential (4>\ c 
classical coupling 
discarded. Notice 



obtained in the previous section ex- 
rmalization countcrterms. Terms pro- 
portional to u 2 or A 2 in the quantum corrections can not 
there are no A or u terms in the clas- 
0). On the other hand, there is a 
g and higher order terms in g can be 
that for m 2 — the result is the same 



as for the charged superfluic 



F e/ (^ C) m = 0)«- 



m 2 m p9 
p 12tt 2 



+ 



8tt 2 



4 



In 



4 



25 



9 At 



(24) 



and therefore there is a fluctuation induced quantum first 
order transition varying m 2 for g ~ u 2 . However, the 
correct results in the present case must be obtained for 
m 2 7^ where superconducting fluctuations are impor- 
tant but not critical. The effective potential to lowest 
order in powers of m 2 , which is also small but finite, is 
given by, 



1 



ef 



2i2 



um 2 



+ 



w 3c 



8tt 2 



^3<: 



hi 



16tH 



0(m 4 ) 



25 
IT 



(25) 



where M 2 is a renormalized magnetic mass parameter 



M 2 = m 2 -^ 
p p 12tt 2 



(26) 



We have now in Eq. (|25J) a new term proportional to um 2 
and since it depends on the cut-off A another counterterm 
is needed. With a minimum counterterm we get, 



V, 



efW3c 



9« 2 (</> 3 



2tt 2 



4! 



u 3c 



8tt 2 



8tt 2 
In 



!J 3c 



In 



- 3 



>3c] 



25 

IT 



(27) 



where other small quadratic terms in the field have been 
included in the renormalized mass M 2 . The coupling to 
the massive superconductivity fluctuations changes dra- 
matically the behavior of the potential since now its sec- 
ond derivative at 4>3 C = is always negative, i.e., the ori- 
gin is always a maximum for any um 2 ^ 0. This coupling 
also gives rise to new minima near <j)^ c = which move 
away from the origin as its strength is increased (Fig.[3J). 
It is interesting that as M 2 is reduced, i.e., as the system 
moves away from the SC-QCP towards the AF-QCP, it 
goes, through a quantum first order transition at M£ 2 , 
from a phase with a reduced value of the antiferromag- 
netic order parameter to another with a larger value of 
</>3c (Figs. and |3J. The AF phase closest to the super- 
conductor is always the small moment one. 

In Fig. the ratio R of the sub-lattice magnetizations 
in the small moment antiferromagnetic (SMAF) and 
large moment antiferromagnetic (LMAF) phases is plot- 
ted as a function of the relevant parameters, namely the 
ratio u 2 /g. Small moment antiferromagnetism (SMAF) 
is a common feature in heavy fermion materials 33 . The 



7 



most common values for the ratio R observed experi- 
mentalljiSi^ are of order 10 -2 , 10~ 3 and correspond to 
values of u 2 / g w 1. We have found that the ratio R is 
independent of the values of g and A being a function 
exclusively of the scaled variable u 2 /g. 

There is an interesting similarity between the phase di- 
agram of the heavy fermion YbRh2Si2, for which a first 
order phase transition from a low moment AF phase to 
a large moment magnetic phase has been observed with 
increasing pressure 3 ^ and the results obtained above. In 
Yb based heavy fermions, pressure (P) acts on the op- 
posite direction it does on Ce compounds decreasing 
the ratio J/W between the Kondo lattice parameters^. 
YbRh 2 Si 2 at P = is a SMAF system with T N = 70mk. 
Pressure increases T/v and at P c w lOGPa there is 
a first order transition to a high moment state with 
fiYb ~ 1.9/iB. On the other hand, negative pressure, 
i. e., expansion of the lattice drives this system to an an- 
tiferromagnetic quantum critical points. It is an excit- 
ing possibility that further expansion of the lattice would 
give rise to superconductivity. 





FIG. 3: New minima appear in the potential for urn 2 7^ 0. 
The effective potential is shown here for two situations: for 



M 2 



Mp , where the first order transition from LMAF to 



SMAF occurs and these two states become degenerate and 
for the spinodal point at which the LMAF becomes unstable 
inside the SMAF phase. 



When the separation A 2 between the classical T = 
critical points is reduced the magnetic moment in the 
SMAF goes continuously to zero. In practice, when the 
separation is small enough, the moment in the SMAF is 
so low that it can be identified with the normal state. 
Then for A 2 sufficiently small, but still large enough to 
avoid the possibility of a direct transition from AF to SC, 
there is a first order AF-normal transition. The spinodal 
points and the energy equivalent to the latent heat can be 
calculated exactly as in the charged superfluid 22,36 since 
the term proportional to um 2 can be neglected. Particu- 
larly, the effective potential at the asymmetric minimum 
close to the first order transition can be written as 



Veffdhc}) « ~M 2 V { 



M 2 



(28) 



u=0 



fit 

AF \f^t 



SC 



M 



u*0 



LMAF 1 SMAF 

"IT 



SC 



0. 



m,p and m 



A 2 is the 
2 are the 



FIG. 4: Phase diagrams for u — and u 7^ 
distance between the AF and SC QCPs 
distances from the point the system is actually probed to the 
AF and SC QCPs respectively. When M 2 w m 2 is reduced, 
i.e., the system moves away from the SC-QCP, there is a first 
order transition at Mp 2 between two AF phases with different 
values of the order parameter. 




FIG. 5: The ratio R = {<f>f M )/(^ M ) between the ground 
state sub-lattice magnetizations in the phases SMAF and 
LMAF as a function of g 2 /u. 



where the critical mass m 2 is 



2 ou / / \2 
*c = 7^ 3 > ■ 



(29) 



Therefore, the latent heat is given by the simple expres- 
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sion 

L h = \m 2 c {fo) 2 . (30) 

A similar phenomenon occurs for the first order SMAF- 
LMAF transition. In this case there is also an associ- 
ated latent heat and the spinodal point corresponding to 
the limit of stability of the large moment phase into the 
SMAF phase is shown in Fig. [31 



VII. QUANTUM BICRITICAL POINT 

In this section we study the possibility of a quantum 
bicritical point (QBP) as shown in Fig The quan- 
tum corrections can be calculated in the paramagnetic 
side and involve terms of higher order in the couplings 
u, A and g. In the limit of small coupling considered in 
this work, these corrections are then of little importance, 
especially for the quantum transition, since the classi- 
cal part becomes much larger. However, mass renormal- 
ization can change the properties of the superconductor 
itself as, for example, it affects the constant k which de- 
termines whether the system is a type I or type Ili£. 

The transition where both phases collapse is continu- 
ous, as obtained classically and the analysis of this dou- 
ble critical point can be done more easily using a quan- 
tum scaling theory^!. Classically this point has SO(n) 
symmetry^, but when taking into account its quantum 
character, as it occurs at T = 0, there are different dy- 
namics related to the distinct phases on both sides of the 
transition. It is easy to convince oneself that in the low 
frequency, long wavelength limit, it is the slow relaxation 
dynamics of the magnetic component with the associated 
z = 2 dynamic exponent that is relevant for the critical 
behavior. For d — 3, assuming that a unique correlation 
length diverges as T is reduced at the QBP, we obtain 
that both the antiferromagnetic (Tjv) and superconduct- 
ing (Ts) critical lines of finite temperature phase transi- 
tions rise with the distance 6 to the QBP, as Tjv,s <x \5\^, 
i.e., with the same shift exponent -0- Furthermore, we 
find ip = vz, where v and z are the correlation length 
and the dynamic exponents associated with the QBP. 
Since the effective dimension rf e // = d + z = 5, we have 
v = 1/2, such that the quantum bicritical crossover ex- 
ponent^ vz — 1. Alternatively, the critical line can be 
written as S(T) = 5(T = 0) + uT 1 /^. For the correlation 
length along the critical trajectory, i.e., for 8(T = 0) = 
scaling yields 

£ ~ \5(T)\- V 5< ^° £,~T- V ^ = T- 1 ' 2 (31) 

so that at the QBP the correlation length diverges with 
decreasing temperature as 1/ VT. The scaling of other 
physical quantities can be predicted similarly, particu- 
larly, the specific heat along this line has a non-Fermi 
liquid, C P (T) cx VT, behavior^ (see Fig.®. 




8=0 

FIG. 6: A quantum bicritical point separating an antiferro- 
magnet from a superconductor. Both lines of finite temper- 
ature phase transitions rise with the same exponent ip. For 
d=3, ij) = vz = 1. 



VIII. CONCLUSIONS 

In this paper we have considered a generalized 
Ginzburg-Landau functional including both supercon- 
ducting and magnetic order parameters to study the mu- 
tual influence of these instabilities in three dimensional 
heavy fermion metals at zero temperature. The quantum 
corrections to this classical functional were obtained us- 
ing the effective potential method up to one loop order. 
This method had to be generalized to take into account 
the dissipative nature of the propagator of the param- 
agnons in the metallic paramagnetic phase close to the 
antiferromagnetic instability. The superconductor prop- 
agator we used corresponds to the simplest generaliza- 
tion of the Ginzburg-Landau classical action to take into 
account dynamic or quantum effects. The calculations 
however can be carried out with different propagators, 
for example, that of Ramazashvili and Coleman? 3 . The 
results will be presented elsewhere. We point out that 
the exact form of the superconductor propagator is ir- 
relevant for the qualitative results of sections VI-A and 
VII. 

Different phase diagrams have been considered accord- 
ing to the distance between the quantum critical points in 
the classical potential. For non-coincident QCPs the sec- 
ond order nature of the quantum transitions obtained at 
the classical level is modified by the quantum corrections, 
changing from continuous to weak first order transitions. 
Hence, the spin fluctuations can change the nature of the 
superconducting transition in the same way that the cou- 
pling to a magnetic field does^S^. For this fluctuation 
induced first order transition we can calculate the equiv- 
alent of the latent heat to transitions at T — and the 
spinodal points. Both indicate that we are dealing with 
a weak first order transition as they are directly related 
to the small coupling constant u. 

We also have found that a metallic antiferromagnetic 
QCP is strongly affected by the presence of supercon- 
ducting fluctuations. These lead to the appearance of a 
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small moment phase which goes through a weak first or- 
der transition to a large moment one as the system moves 
away from the superconductor instability. 

Finally, in the case of coincident quantum critical 
points the transition remains second order and the finite 
temperature behavior can be extracted from the scaling 
properties of the quantum bicritical point. 

The effects of the competition between the AF and SC 
instabilities are highly non-trivial since a positive cou- 
pling which classically tends to avoid the coexistence of 
these states leads when quantum corrections are consid- 
ered to symmetry breaking in the normal state and in- 
creases the region of the ordered phases (SC and AF) in 
the phase diagram. 

In recent years increasing experimental evidence is be- 
ing gathered that points to existence of intrinsic inhomo- 
geneities close to the T = antiferromagnetic instabil- 



ity in heavy fermionsMi. Coexistence of small ordered 
moments with paramagnetic regions, coexistence be- 
tween superconductivity and paramagnetism have been 
observed. The present approach provides an appealing 
model for this type of phenomena. They arise in our the- 
ory as a direct consequence of the competition between 
different order parameters which give rise, through quan- 
tum corrections, to weak first order transitions with the 
associated physics of metastability, spinodal points and 
phase coexistence. 
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